Preprint typeset in JHBP style - PAPER VERSION 



*** 



Density-metric unimodular gravity: vacuum maximal 
symmetry 



Amir H. Abbassi 

Department of Physics, School of Sciences, Tarbiat Modares University, 
P.O.Box 14155-4838, Tehran, Iran. 
i?-maiZ;ahabbasi@modares . ac . ir 

Amir M. Abbassi 

Department of Physics, University of Tehran, 
P.O.Box 14155-6455, Tehran, Iran. 
i?-maiZ;amabasi@khayam . ut . ac . ir 



Abstract: We have investigated the vacuum maximally symmetric solutions of recently 
proposed density-metric unimodular gravity theory, the results are widely different from 
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1. Introduction 

Unimodular gravity is a very modest modification of Einstein general theory of relativity 
which leads to the same field equations with the cosmological constant emerging as an 
integration constant. In unimodular theory one apparently reduces the dynamical compo- 
nents of the metric by one. This is accomplished by discharging the determinant of the 
metric as dynamical variable. Conventionally this demand is granted by imposing = 1 
as a constraint [1-9]. However, it was shown that one indeed recovers the same number 
of dynamical degrees of freedom as in GR via a tertiary constraint [10]. Even though a 
covariant expression on variation of the metric, g^v^g^v = emerges from the unitary 
condition \g\ = 1, but the primary relation itself is not established on a covariant form. 
Determinant of the metric is a scalar density and cannot be kept equal to one in each 
arbitrary coordinate system. Evidently, this is not compatible with our notion of general 
covariance. This disadvantage has been removed in density-metric unimodular gravity pro- 
posed recently[ll]. In metric theories of gravity the gravitational field is represented by 
the spacetime metric which is a symmetric tensor field of rank (0, 2) and signature of +2 
defined on a four dimensional manifold. In density-metric theory the gravitational field 
is represented by a symmetric tensor density field of rank (0,2) and weight —\ which its 
determinant is equal to one by definition. Since the determinant of a tensor density of 
weight — ^ in a four dimensional spacetime is a scalar under coordinate transformations, so 
the unitary condition may be imposed consistently with the principle of general covariance. 
Using this density metric we may equip the manifold by a specific connection. It turns out 
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that introducing torsion-free and metric compatibility conditions are not sufficient to fix 
uniquely all components of the connection. Actually, the components of trace of connection 
remain undetermined. Hence the suppositions are: 

1 , , 

9(Ui/ — 9i^/u Symmetric tensor density of weight — — (l-l) 
|g| = 1 Unitary condition (1-2) 
r^. = re^ Torsion free (1.3) 
^P^ni/ ~ Metric compatible (1-4) 

By expanding out the equation (1.4) for three different permutations of indices and sub- 
tracting the second and third of these form the first we obtain: 

^\^'p^^I^^'pA-V\.) (1-5) 
Contraction of indices A and in (1.5) gives 

ri = ^9^''5,9;,, + r^, (1.6) 

Obviously (1.6) brings us to 

9^''^.3ap = (1.7) 

So four relations of (1.7) reduce ten independent components of density-metric to six. 
Therefore components of connection typically are functions of six independent components 
of density- metric and four independent components of the trace of connection Fp^. Us- 
ing the connection(1.5), we may introduce Riemann tensors and Ricci tensor and scalar 
respectively by: 

= d^'^U - d^^Kc + Kx'^la - r^A^^. (1-8) 

Rau = (1.10) 
R = fRau (1.11) 

Here (1.8) is a tensor of rank (1,3), (1.9) is a tensor density of rank (0,4) and weight — ^, 
(1.10) is a tensor of rank (0,2) and (1.11) is a scalar density of weight Let us notice 
some algebraic properties of Riemann tensor. It is antisymmetric in its first two indices 

(1.12) 

and it is not invariant under interchange of the first pair of indices with the second: 

Rpap,v 7^ Rjxvpa (-'-•-'-^) 

The sum of the cyclic permutations of the last three indices of Riemann tensor vanishes: 

Rpafxv ~\~ Rpixvcr ~l" Rpuap, ~ (1.14) 
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The Bianchi identity holds too. The sum of cychc permutation of the first three indices of 
covariant derivatives of the Riemann tensor is zero: 



Rpau;X + Rpa\n;v + RpauX;n = (1-15) 

The Ricci tensor associated with the Christoffel connection is symmetric but in our case 
this is not generally true. Twice contraction of Bianchi identity (1.15) with g^'* and q'^'^ 
plus doing some calculations in Riemann normal coordinates result in: 

V^l^i^A. + iRaX - l^axR) = (1-16) 

In density-metric unimodular gravity figured that was the simplest possible choice for 
a Lagrangian density and proposed 



= J KR^dx"^ (1.17) 



for the action of gravitational fields in the absence of matter and other kinds of sources. 
Variation of the action (1.17) with respect to F^^ leads to the field equation 

VAii; = (1.18) 

Also variations of the action (1.17) with respect to g'^'^ and applying the unimodular con- 
dition 

by the method of Lagrange undetermined multipliers and inserting field equation (1.18) 
lead to 

%H-iV^ = (1-20) 

We should mention that density-metric is not eligible to carry out the whole tasks of a 
common metric. First its operation on two vectors gives a scalar density instead of a 
scalar. Second the density-metric has six independent components. In accordance with 
affine theories of gravity we may equip the spacetime with a metric by defining 

9t^u = R%^ (1-21) 

where R is given by eq.(l.ll) and of course should be nonzero. In addition it depends 
on four independent components of trace of connection,Fpju. Thus g^i, depends on ten 
independent dynamical variables. Given the equation (1.21) in principle, we can work out 
density-metric and trace of connection in terms of ten independent components of g^i, and 
formally write: 

%u = %Agap) , r;;^ = r^^(5a/3) (1.22) 

We have already worked out the solution of this theory for vacuum spherically symmetric 
space-time: 

ds"^ = B{r)dt^ - A{r)dr'^ - r^{de^ + sir?ed4>'^) (1.23) 
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with 



and 



Bir) = (1 + (1 + + Kr\l + -^)-^] (1.24) 

..X , C\ 2GM-AC' C\ 1 . 2. 1. 1 . X 

A{r) = (1 + ^ (1 + + Ar^(l + ^)-.]-^ (1.25) 

In these relations M is the mass of point object which is the source of sphericahy symmetry 
and A has dimensions [L~^] where may be identified as cosmological constant. C is a new 
constant which appears in this theory with dimensions [L~^]. 

The constant C has a physical interpretation, a convenient way to show its real nature is 
when we manipulate the post-post-Newtonian formalism ; just as the same way that the 
Schwarzschild radius is fixed by comparison with Newtonian limits. It should be noticed 
that a spherically symmetric matter distribution has a total mass, but certainly has a 
vanishing quadrupole moment. Such a matter distribution could have three equal principal 
moments of inertia I. It seems very probable that C" oc ^ . In this way C is an individual 
property of the source object. Observational data from double pulsar systems can be used 
to verify this issue. Here we can offer an upper bound for it. 

In the next section we explicitly work out these results and correctly predict what are 
observed in three classical tests of gravitation. 

2. Classical Tests 

For this purpose we may choose A = 0. The geodesic equations imply that the quantity 

is constant along the path. For massive particles we typically choose A so that e = 1 . While 
for massless particles we always have e = 0. Invariance under time translations leads to 
conservation of energy while invariance under spatial rotations leads to conservation of the 
three components of angular momentum. Conservation of direction of angular momentum 
means that particle moves on a plane. So we may choose 9 = ^. The two remaining Killing 
vectors correspond to the energy and magnitude of angular momentum. The energy arises 
from time-like Killing vector {dtY = (1)0,0,0) while the Killing vector whose conserved 
quantity is the magnitude of angular momentum is (5,/,)'* = (0, 0, 0, 1). The two conserved 
quantities are 

C\ 2,_ 2GM , C\ i,dt 
and 

-^^f (.3, 
Prom eq.(2.1)by using eqs.(2.2) and (2.3) we can arrive at the following result: 
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Then expanding in powers of r, we shall content ourselves here with only the lowest order 
of approximation in eq.(2.4). So we have 

+ + e)-E\l + -—) + ^ = (2.5) 

Now we are prepared to discuss the classical tests one by one. 

1) Precession of Perihelia: With e = 1 in eq.(2.5) and multiplying it by (^x)"^ ~ 
defining x = -^j^ we have 

_ 2. + - — + -^(2C' - — )x3 = _ - _ (2.6) 
Differentiating eq.(2.6) with respect to gives 

= + (2.7, 

After some manipulation [12] we obtain that the perihelion advances by an angle 

On the other hand we have 

= 1 + — po— (e^ - 1) and ^ GM(1 - e'^)a (2.9) 

where e is the eccentricity and a is the semi-major axis of the ellipse. Putting eq.(2.9) in 
eq.(2.8) yields 

GttGM C AC 

{1- e2)a^ ^ 3G2M2(1 - e'^)a ^ 3GM(1 - e^)a'^' ^ ' 

The third term in the parenthesis is smaller than the second term by a factor and may 
be neglected. 

'^^=(r3^" + 3G^M^(l-e> ) 
Observational data put an upper bound on the C. 

2) Deflection of Light: Putting e = in eq.(2.5) and multiplying it by (^)~^ = ^ and 
defining u = ^ we have 

rill F2 

(^)2 + „2_2GfMu3 = f_(i + 8CV) (2.12) 
dcp 

Differentiation of eq.(2.12) with respect to (f> makes 

5+- = (3GM + i^).^ (2.13) 
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By some straight forward calculations [13] the deflection angle is 



4 4C' E'^ 

where ro is the distance of closest approach to the origin. Since j; is ^ then eq.(2.14) 
becomes 

^ 4GM 16C" 

S = h — ^ (2.15) 

ro ro^ 

An upper bound on C may be obtained by comparison with observational data. 

3) Gravitational Redshift: Let us Consider an observer with four- velocity who is station- 
ary = 0), we have: 

U° = (B(r))-i (2.16) 

Any such observer measures the frequency of a photon following along a null geodesic x^i'^) 
to be 

dr^ 1 dt 

„ = = (2.17) 

By considering cq.(2.2) and (2.2) we conclude that for a photon emitted at ri and observed 
at r2 the observed and emitted frequencies are related by 



In the interval A~2 > r » 2GM and C s we get 

^ = 1 - GM(1 - 1) - ^(i, - i,) + |(r? - 4) (2.19) 
This is also consistent with observational constraints. 



Now we want to work out the solutions of these field equations for the vacuum maximal 
symmetry space-times. In other words our goal is to illustrate the analogue of Minkowski , 
deSitter and anti- deSitter space-times in the density metric unimodular theory of gravity. 
The deSitter-Schwarzschild analogue in static form may be obtained by putting M = C = 
in eqs.(1.24)and (1.25). 



3. McLximal Symmetry 

A rigorous mathematical analysis may be applied by taking N dimensional manifold 
equipped with a density metric and a trace of connection. Assuming the density metric of 
rank (0,2) and weight — ^ being form invariant under a given coordinate transformation 
generated by a Killing vector ^'^ leads to: 

V^e. + V.^^ - ^S^^VaC^ = (3.1) 
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This means that ^ f_iCiy~'zi^i^^ is an anti-symmetric tensor. Since the trace of connection 
is taken as dynamical variable the condition to be form invariant provide us: 

v.VaC^ = ^'(^^r^A - - ^^A^' (3.2) 



It can be easily shown that 



V.VpCm = -Rapi^^X - ^9api?M?^ + \^aA>^^' + ^V^-A^^ (3.3) 

Following the same argument as given by Weinberg[14], we may conclude that the manifold 
should have ^^^^"^-^ Killing vectors to be maximally symmetric. We start from general 
formula for commutator of covariant derivatives of tensors 

V.V.(V^Cp - - v.v,(v^ep - l^^p^vC) = 

-RpaA^^^^x - i9A^V,e'') - <.(Va?p - \ixpV,C) (3.4) 



Then apply eqs.(3.1),(3.2) and (3.3) into eq.(3.4). It can be shown that for a maximally 
symmetric manifold we should have: 

Ri,u = (3.5) 

This is an important result, since it causes the Riemann curvature tensor to resume the 
same algebraic property just as Riemann geometry, i.e. if Christoffel connections were 
used. Following Weinberg's procedure we obtain 

{N - l)Rxpau = Rup%Xa " Rap^Xu (3.6) 

and ^ 

Rap = J^^apR (3.7) 

Putting (3.7) into (1.16) gives 

VxR = (3.8) 
Of course this does not mean R is constant since it is a density scalar. 

4. Vacuum Cosmic Dynamics 

The homogeneity and isotropy of a space-time imply that it is maximally symmetric. In 
CO- moving coordinates the metrics are given by Robertson- Walker line element, 

ds^ =dt^-a\t)!^j^^+r^{de^ + sm^edcl>^)^ (4.1) 

To prepare the appropriate density metric consider the following scenario: Multiply the 
components of Robertson- Walker metric by inverse of-^/j^ , where|g'| is determinant of the 
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metric. This results in a density tensor of rank (0,2) with the desired weig ht -i . The 
nonzero components of the new density metric are, 

_ (1 — A;r^)3 _ a5 

3 1 ' QfT" 1 3~ 

a^rsin^O r sina ^(1 — fcr^)! 

aarfl — A;r^)4 i . i , 9si x 

%ee = '—,^^ = a2rsm2e{l-kr^)4 (4.2) 

sin 2 

and the nonzero components of its inverse are, 

aarsinz^ r sin 2 9 {1 — kr'^) 4, 

9 = ~~. I7T > 9 



I — kr'^)^ a2 

0"' = ■ f'!.,. . = ■ , ' , (4.3) 

Q,2r(l — A;r^)4 02 r sina 6'(1 — A;r^) 4 

In terms of the unknown components of trace of connection, the nonzero components of 
the connection are determined by using eqs.(4.2),(4.3) and (1.5). The results are presented 
in A.l. The off diagonal components of the field equation (1.10) simply are: 

R^^. = for n^u (4.4) 

Inserting the connection from A.l into eqs.(1.8), and (1.10), then Rj-t = gives 

Equation (4.5) derives a manifestly unique expression for the radial component of trace of 
connections. 

r^' = - + (4.6) 

Similarly the polar component of trace of connection will be obtained by considering the 
equation Rtg = 0. It gives : 

(r^* + ^)(r^.-cot0)-2|r^, = o (4.7) 

Eq. (4.7) imposes a unique expression for polar component of connection trace 

^% = cote (4.8) 

The azimTith component of trace of connection will emerge immediately by considering 
Rt(l> = . This simply reads 

+ ^) = (4.9) 
Therefore eq.(4.9) at last determines azimuth component of connection trace 



K<P = (4.10) 
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The other off diagonal components of field equations, Rj-g = R.^.^ = Rg^f, = , are 
satisfied automatically by the given equations(4.6),(4.8)and(4.10). Plugging equations 
(4.6),(4.8)and (4.10) into the given components of the connection in A.l reduce the non- 
vanishing components of the connection. They are given in A. 2. These results may be used 
to compute the diagonal components of Ricci tensor: 

^« = i(^ + ^^* + ^^*^ ^'-''^ 

Res = -2kr' - iaV^P^ + + 2r% + 4^r^, + T%'] (4.13) 
R^^ = R0esiT?e (4.14) 

Let us write the diagonal components of the field equation (1.20) : 

^3"i?tt - l^'^Rrr - i''Roe = (4.15) 

l^'^Rrr - y^Rtt - %''Ree = (4.16) 

^''Ree - l^'^Rtt - l^'^Rrr = (4.17) 

The azimuth component of eq.(1.20) reproduces the same equation as eq.(4.17). To achieve 
the final result it is sufficient to substitute eqs. (4.11)-(4.13) into eqs. (4.15)-(4.17) which 
all merely get to the same equation as, 

The field equation (1.18) does not express a new equation and indeed relation (4.18) already 
satisfies equation (1.18). In the case of Minkowski space-time, a ,is constant and k is equal 
to zero. These are in agreement with relation (4.18) if we distinguish, 

r^t = (4.19) 

According to arguments expounded in [11], concerning gauge fixing, we regard that generally 
relation (4.19) rigorously holds. This is consistent with our result for vacuum spherical 
symmetry and M = presented in the mentioned reference. Finally substituting relation 
(4.19) into eq.(4.18) yields 

a l.d.n 4A; „ x 

- - 7 - - ^ = 4.20 
a 4 a 

It is relevant to compare the obtained equation (4.20) with its counterpart in general theory 
of relativity which is Priedmann equation. We have [15], 

d2 = ^a^ -k (4.21) 

2^ + (^)2 + A-A = (4.22) 
a a 
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Eliminating A between equations (4.21) and(4.22) we get, 



- i-r -4 = (4-23) 



a a 

It is interesting to notice that in eqs. (4.20) and (4.23) identical terms appear with merely 
different numerical coefficients and both equations are invariant under time reversal trans- 
formation. So it is plausible while they have probably the same roots we predict to obtain 
results with different behavior than the existing models. We are going to solve the eq.(4.20). 
To do this taskjwe assume in accordance with eq.(4.21) it has a solution of general form 

d2 = /(a) (4.24) 

Differentiating eq.(4.24) with respect to time and taking a ^ yields 

« = 5^ (4.25) 
2 da ^ ' 

Plugging eqs. (4.24) and (4.25) into eq.(4.20) gives 
Eq.(4.24) has the general solution, 

= /(a) = AVa - 16A; (4.27) 

where A is constant of integration and we have intended to call it so to show its resemblance 
with cosmological constant in Priedmann model. Eq.(4.27) may be integrated and has the 
general solution 

±t= -"^^^t"' VX^^ + C (4.28) 

where C is the second constant of integration and may be put equal to zero without loose 
of generality. To find scale factor as a function of time from eq.(4.28) we will face with a 
cubic equation. In the next section we will discuss some features of the solutions of this 
equation in different situations. 

5. Empty Cosmological Models 

Reality condition imposed on (4.27) and (4.28) imply the following limits on the domains 
of a and t without solving the equations. 

16 

a > {^f t G R deSitter (5.1) 

a > teR deSitter (5.2) 

2 16 

a > |t| > -(-7-)^ deSitter (5.3) 

no solution (5-4) 

a = const t & R Minkowski (5.5) 



I 


A > 


k 


= 1 


II 


A > 


k 


= 


III 


A > 


k 


= -1 


IV 


A = 


k 


= 1 


V 


A = 


k 


= 
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VI 


A = 


k = 


-1 


VII 


A < 


k = 


+1 


VIII 


A < 


k = 





IX 


A < 


k = 


-1 



aoct teR Milne (5.6) 

no solution (5.7) 

no solution (5.8) 

16 2 16 

< a < (-r)^ 1*1 < -(-r)^ anti-deSitter (5.9) 

Now we express some points for each case in more detail. 

I. A > 0,A; = 1 

All closed models with positive lambda have a lower bound on a equal to o^m = 
To find a as a function of time from eq.(4.28) we face with a cubic equation which its 
discriminant is positive for this model. This means it merely has one real root, it is: 



a = {^-lm + 9( A)4i2 + + (1 + 9( A)4t2)2 



+ \J 1 + H^rt' - ^-1 + (1 + 9(^)4^2^^ - 1]}2 (5.10) 



Models of class I possess event horizon which can be find by considering 

r« dr _ /•+~ _dt 
Jo VI -r^ Jto «( 



(5.11) 



We infer from eq.(5.10) that a +oo as t ^ +oo. Taking a{to) = oq, eq.(5.11) gives 

sin-\rH)= r°°— (5.12) 

Jto «« 

Calculating the integral in eq.(5.12) by inserting eq.(4.27), we will have 



rn = cos{tan-^ y - 1) (5.13) 



II. A > 0, A; = 

In flat models with positive lambda eq.(4.27) reduces to 



^ = ±VA (5.14) 



and integration of eq.(5.14) with respect to time gives 

a = {\^\t\)t (5.15) 

The integration constant has been chosen so that wc have a, = at t = 0. While 

Hubble and deceleration parameters in Fricdmann empty flat model are constant, 
eq.(5.15) implies 

H = ^ = ^ (5.16) 

a St ^ ^ 
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and 

Deceleration parameter eq.(5.17) stands for an accelerating Universe which its accel- 
eration rate decreases by passing of time. 

The exponential growth of deSitter space in models / and // is replaced by a power 
law which practically is more plausible to match the physical reality. This behavior 
is in no way inflationary. Model II also possesses an event horizon 

f^°^ dt f^^ 1 ,^0^4 , 1 
rH= -77S= -{T^'dt = — 5.18 
J to J to «o t 3ao 

where qq = a{to). It should be noted this model is singular at t = 0. 

III. A > 0,fc = -1 

Eq.(4.28) for an open model with positive lambda has a real root 



+ ^-1 + H^rt' - y -1 + (-1 + H^m' - i^-m 

which is defined merely for \t\ > While it is singular at t = 

defined for long intervals of time. For this reason this model can not present a 

physically acceptable space-time. 

IV. A = 0, A; = 1 

There is no solution for a closed model with zero lambda. 

V. A = 0,A; = 

The flat model with zero lambda presents the well-known Minkowski space. 

VI. A = 0,A; = -1 

The open model with zero lambda presents Milne space. 

VII. A < 0, A; = +1 

There is no solution for closed models with negative lambda. 

VIII. A < 0, A; = 

There is no solution for flat models with negative lambda. 

IX. A < 0,A; = -1 

In this case eq.(4.28) has a negative lambda has a negative discriminant which means 
there are three real roots. Since this does not give a unique solution and the definition 
of function is not fulfilled, so should be rejected. 
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6. Concluding Remarks 



1- There exists no solution with oscillating nature. 

2- Models // and /// are singular and the nature of their singularities could not be 
explored by constructing a scalar of Riemann tensor. It is illustrative to compute the 
density scalar R by applying eqs.(4.25) and (4.26), we get 

„ 3 r sin2 ^ 

R = -- tA 6.1 

4(1-A;r2)4 

The action may be rewritten by inserting eq.(6.1) into eq.(1.17) 

This does reminisce a non-fiat geometry with constant curvature. So we may hope 
with some certain that singularity is non-intrinsic. 

3- The Riemann tensor made of the metric (4.2) while satisfying (4.27) is consistent with 
eq.(3.6). So we may say with confidence that the obtained space-times are maximally 
symmetric. 

4- Open models admit merely A = 0. 

5- The inflationary behavior of dcSittcr models are replaced by a power law character 
which reserve to be physically interpreted. 

6- The event horizon rn given by eq.(5.18)goes to infinity if to approaches to zero. We 
may hope while we have improved the inflationary scenario also have resolved the 
existence of event horizon. 

A. 1 

Non-zero components of the connection are: 

r?. = ^(r^.-3^) , ^\r = \{T^,r-\--^) , r*, = i(r^,-cot^) , r^, = ir^, 

■pt frp ^\ -pt '^^^^ /pP ^\ T* a^r^sin^^ p d 

^ 1 - fer^ _ ^ = -(rp -\--) r*" = -(rp --4- ) 

tt AO K'- or „) /I „2 ' Pt^ „J ^ '-rr a\ pr „ ^ -, i„„2 > 



Ao? ^ P^ r' Ao? ' A^ P^^ a' ' A^ p"" r^l-kr"^ 

Vie = \{rp^- cote) , K^ = \^U ' ne = -lAi-kr'){rp, + ^) + \kr^ , r^^^ = sm^ eri^ 

^tt = -^,iKe-^-'(^) ' ^^* = i^^P* + ^ ' rf^ = - 4r2(i-fc,2) (r^.-^otg) 
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P<?!> _ p<f> p<i!> _ _(-pP r<t> = £^ r*^ 

** 4a2r2sin2 0' 4^*^^^' r-r 4r2(l-fcr2)sin2e ' 't^ 

1 1 

= ^ , rf. = -(r^. + Scot^) , , = -r^,. 

4 sin 9 't'" 4 ^ P^ 4 p<p 

(•) represents differentiation with respect to time. 

A. 2 

The non-zero components of the connection after inserting eqs.(4.6),(4.8) and (4.10) in the 
given components of connection in A.l are: 

= \{r% - A , = + 1) , = \a\Hr% + ^) , r^, = r*, sin2 e 

^et = lK + l) ' ^er = l ' r^^^ = -sin^cos0 
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